SUMMATION OF HYPERHARMONIC SERIES 



ISTVAN MEZO 



Abstract. We shall show that the sum of the series formed by 
the so-called hyperharmonic numbers can be expressed in terms 
of the Riemann zeta function. More exactly, we give summation 
formula for the general hyperharmonic series. 



1. Hyperharmonic numbers 

Introduction. In 1996, J. H. Conway and R. K. Guy in |CG] have 
defined the notion of hyperharmonic numbers. 

The n-th harmonic number is the n-th partial sum of the harmonic 
series: 

n 1 

k=l 

rW ._ 



H\ := H n , and for all r > 1 let 



(r-l) 
k 

k=l 



be the n-th hyperharmonic number of order r. These numbers can be 
expressed by binomial coefficients and ordinary harmonic numbers: 

a) Fr = (^l~ 1 )(^ + ,-i-^,- 1 ). 

It turned out that the hyperharmonic numbers have many combi- 
natorial connections. To present this fact, we need to introduce the 
notion of r-Stirling numbers. Getting deeper insight, see [BGG] and 
the references given there. 

r-Stirling numbers. [™] r is the number of permutations of the set 
{1, . . . , n} having k disjoint, non-empty cycles, in which the elements 
1 through r are restricted to appear in different cycles. 

The following identity integrates the hyperharmonic- and the r- 
Stirling numbers. 



rn+rl 
Lr+lJ , 



H {r) . 



nl 
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This equality will be used in the special case r = 1 |GKPj : 

rn+ll Tn+ll 

(2) ^ : = ^ = H »- 

2. Results up to the present 
Our goal is to determine the sum of the series has the form 

oo n ( r ) 

71=1 

for all r > 2 and possible m. 

To do this, we should determine the asymptotic behaviour of hyper- 
harmonic numbers. In the paper [M] there is a skimped approximation 
which helps us to get convergence theorems for hyperharmonic series. 
Namely, we have that 

1 „r- 1 „ jj(r) „ ^ (2r) r _ r 



:7l' < < 



r! n 2 (r-1)! ' 

for all n G N and r > 2. According to this result, the followings were 
proved: 

oo_ „(r) 
-l-ln 



77.' 
n=l 



C(s + 1) A/ft r) ,3 (2rf 
^T- < ^^ < 2(73 1 )T C(S) (S>1) ' 

n=l v ; 



oo __1 

n 

oo, 



n=l 

9 (V — 1 V r> r ~ s 

- , , ^ )<y^T<r!((^l) (s>2), 
3 f2r) r ^ f/( r ') 

where C is the Riemann zeta function and r > 2. 

3. Asymptotic approximation 

To have the exact asymptotic behaviour of hyperharmonic numbers 
we need the following inequality from [CGj . 

(3) 2 (w + l) + ln 0)+7<#n< — + ln(n) + 7 (n G N), 

where 7 = 0.5772 ... is the Euler-Mascheroni constant. 
Now we formulate our first result. 



Lemma 1. For all n GN and for a fixed order r > 2 we have 

1 

(r-1)! 
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that is, the quotient of the left and right hand side tends to 1. 
Proof. The binomial coefficient in ([1]) has the form 

'n + r — 1\ (n + r — l)\ 1 , , , , 

1 - (n + l)(n + 2) -• -(n + r- 1). 



r — 1 J (r — l)\n\ (r — 1)! 

It means that for a fixed order r 

n + r - 1\ 1 -, 

~ t rr™ • 

r-1 / (r-1)! 

For the convenience let us introduce the abbreviation t := r — 1. We 
should estimate the factor if n+ t — H t in (CD). According to ([3]), we get 
that 

H "» - < + ln( " + f) + 7 - 2(tTT) - ln(i) - ''■ 

Since 2 , 1 + ^ < 2 (i+t) an< ^ m ( n + — ^ n (t) = m (*^)> so 

H n+t -H t <ln(n + t). 
The lower estimation can be deduced as follows 

- g ' > 2(n + t + 1) + + '> + 7 - Tt - '""> - 7 > 

> ln(n + t) - ln(i) - - = ln(n + t)- ki(ty/e). 
From these we get that 

ln(n + 1) - ln(iVe) < - H t < ln(n + *), 

whence _ 

1 ffi(Ve) ^ # w+f - fl, < 1 

ln(n + t) ln(n + t) 
The limit of the left-hand side formula is 1 as n tends to infinity. There- 
fore (remember that t = r — 1) 

H n+r -i - H r -x ~ ln(n + r - 1) ~ ln(ra). 

Collecting the results above we get the statement of the Lemma. □ 

Corollary 2. The following series are convergent 

co (r) 

< +oo, 

n m 

n=l 



whenever m > r. 

Proof. Because of Lemma 1, 



y ^ < c y nT ~ x Hn) < c f ^ < +oo 

n=l n=l n=l 



under the assumption m > r + 1. □ 
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4. Generating functions, Euler sums and Hypergeometric 

Series 

In this section we introduce the notions needed in the proof. 
Generating functions. Let (a n ) nS N be a real sequence. Then the 
function 

oo 

a n z' 

n=0 

is called the generating function of (a n ) n€ ^. If a n = H n we get that 

(see [DEB USD]) 

*— ' 1 — z 

n=0 

and in general 



/(*) : =E 



(4) e#<w 



ln(l - z) 

' (1 - zY 

n=0 v 7 



r) n 



The generating function 

(5) -L(-ln(l-«)) m = E " 



771! 

n=l 



n! 



can be found in [UKPl lB]. 

The well known polylogarithm functions can also be considered as 
generating functions belong to a n = \ (for a fixed k). 



oo n 



n=l 



The last remarkable function needed by us is 



1 — 2;) fe V 

y n=0 



Euler sums. The general Euler sum is an infinite sum whose general 
term is a product of harmonic numbers divided by some power of n, 
see the comprehensive paper [FSj . The sum 



m— 2 



E ^ = 2 (m + 2)c(m + x ) - E c(™ - *oc(* + 1) 

n=l fc=l 

was derived by Euler (see [BB] and the references given there). Related 
series were studied by De Doelder in [dDj and Shen [S], for instance. 

Hypergeometric series. The Pochhammer symbol is defined by 
the formula 

(6) (x) n = x(x + 1) • • • (x + n — 1), 
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with special cases (l) n = n\ and (x)i = x. The definition of the 
hypergeometric function (or hypergeometric series) is the following: 



TP 
n A m 



0>l, 0,2, • • • , CL n 

bi, b 2 , . . . , b m 



E 

k=0 



(ai)fe(fl2)fe ■ • • (Ofe Z k 

{W)k{h)k ■ ■ ■ {b m ) k k\ ' 



This function will appear in the sum of the hyperharmonic numbers. 
We shall need one more statement. 



Lemma 3. We have 



(1 — z)z z 



and for all 2 < r e N 

J {l-z)z r 
or, equivalently, 

ln(z 



\n{z) 



-dz 



hi(z) 



1 



(1 - z)z r ~ 1 (r - l)^- 1 (r - l) 2 z 



k=l v 7 



(1 - z)^ 

up to additive constants 

Proof. The definition of Li 2 readily gives that 

\n(z) 



Li' 2 (l-z) = T 



\n 2 {z) 



' Hz) 



Moreover, 

whence 

Li' (1 - z) + 

The first statement is proved. The second one also can be deduced by 
differentiation. The derivative of the right-hand side has the form 

\n(z) (r -l)z r ~ 2 - (r - l) 2 \n(z)z r ~ 2 -(r - 1) 



2 v ; 



z\n(z) + (1 - z)\n(z) \n(z) 



(l-z)z r - 1 (r-l) 2 ^- 1 ) 2 

zln(2) z r — (r — 1) ln(z)z' 



(r - 1) V 



1 

+ 



(l-z)^ (r-l)(z r ) 2 (r-l)z r 

z r+1 \n(z)(r - 1) - z r (l -z) + (l-z)(r- 1) ln(z),z r + z r (l - z) 

(r-l)z 2r (l-z) 
z\n(z) + (1 - z)\n(z) \n(z) 



z r (l - z) 



Z^l-Z) 1 



as we want. 



□ 
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5. The summation formula 
For the sake of simplicity, we introduce the notations 



S(r,m) := 



oo „( r ) 
11 n 



n=l 



rv 



and 



B(k, m) 



1, 1, 

2, 2, 



., 1, fc + 1 
•, 2 

After these introductory steps we are ready to prove the main theo- 
rem. 



Theorem 4. If r > 2 and m > r + 1, then 



r-i l 

S(r, m) = S(l, m) + — [S(k, m — 1) — B(k, m)) . 



k=l 



Proof. We begin with the generating function (j3J). Division with z and 
integration gives that 



oo rr( r ) 
J-J-n 



E 

n=l 



n 



An integral transformation gives that 



ln(l - z) 

z(l - zY 



dz 



ln(l - z) 
z(l- zY 



ln(z) 



'l-z)z r 



dz. 



-dz 



fc=i v v 

According to (jlj) and (JSJ) one can write 

n 



ln(l - 2) 1 



z) k k 2 {l-z) k J ' 



n=l n=l 
r— 1 / 1 00 00 * 

E (-»E^4E 

fe=l \ n=0 n=0 v 



Z 

n ! 



n + k — 1 

n 



It is obvious that this series is divergent but the generating functions 
works well without any restriction. 

Let us deal with the second term. The Stirling numbers satisfy the 
recurrence relation 



n 


= (n-l) 


n — 1 




n 


- 1" 


A 


k 


+ 


_k 






- 1 



From this 



~n + r 






+ 


n 




n 




= n 






= n 




2 




2_ 


1 




_2_ 



(n > 0). 



+ (n-l)! (n>0). 



Now, 
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can be rewritten as follows 



1 

+ -. 

n 



1 


~n + 1 


1 


n 


n! 


2 


(n-1)! 


_2_ 



Division with n and rearrangement give that 



1 

n\ 

Therefore the second sum is 

oo 



H„ 



n 7r 



E 

n=l 



Z 

n\ 



yHn z r, 
n=l 



n=l 



Since the last member equals to Li 2 (z), it cancels the first member of 
the sum above. Hence 



oo n (r) 
-l-ln 



E 



n=l 
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OO jj J — 1 / oo 1 CO 

n=l fe=l \ n=0 n=0 



n + A; — 1 



An easy induction shows that (after dividing with z, integrating, and 
repeating these steps (m — l)-times and finally substituting z = 1) 
(7) 

1 y (n k - I \ I 

^2 



oo (r) 
tin 



rv 



5(1, m) 



E fr< 

fc=i \ 



r,m — 1) 



E 

n=l fc=l \ n=l 

The last step is the transformation of the last member. 



n 



n 



m—l 



E 



n=l 



n + k — 1 

n 



n 



1 oo 
n-l L^i 



(n + k-l)\ 1 
^ n\(k - 1)1 n m - 1 

n=l v ' 



1 oo 1 

— V (n + + 2) • • - (n + k - 1) r 

{k-l)\^ 'n m - 1 



1 x 



(k — 1)! ^ n" 

because of the definition of the Pochhammer symbol in formula 
On the other hand, the definition of B(k,m) yields that 

(n!) m (* + !)„ ^ 1 (fc + l)„ 



n + l) 



rv. 



n=0 y ' n=0 

The next conversion should be applied: 
fc!(fc+l) n (k + n)\ 



t — 

^ (n + l 



rv. 



rv. rv. 
It means that the equality 



(n + l)(n + 2) ■ • • (n + fc) = (n + l) fc . 



(8) 
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holds. That is, 

n=l v 7 

Considering this and (J7J) the result follows. 

6. Tabular of the low-order sums 



□ 



In the following tabulars we collect the low-order results of the Sum- 
mation Theorem. We used the following identities which can be easily 
derived from (ED and (ED). 



B(l,m) = ({m- 1) 
B{2, m) 

B(3,m] 



- (C(r^- 1) 



C(m-2)) 



A computation with the mathematical package Maple shows that 
an improved accuracy can be reached using Theorem 1, in spite of 
calculating the series term-by-term. Namely, the sum 



10° 

E 

n=l 



H. 



1.310972037, 



TV 



and takes about 500 seconds on an average personal computer, while 
the "closed form" shown below gives immediately the closer value 



oo (4) 



E 

n=l 



n" 



TV 

540 



7T 

810 



C(3) - ^C(3) 



with ten digits accuracy. 



216 " v y 36 
» 1.310990854, 

5(2, m) 



^C(3) 2 + yC(5) 



Power of n 


Closed form 


Approx. value 


m = 3 


£-£ + 2C(3) 


2.112083781 


m = 4 


f^ + 3C(5)-C(3)(l + 4) 


1.284326055 


m = 5 


^-^-^C(3) 2 + 3C(5)-^C(3) 


1.109035642 


m = 6 


£ + 4C(7) - gC(3) - K(3) 2 ^ 
C(5) (l + f ) 


1.047657410 


m = 7 


& - m ~ C(5)C(3) + 4C(7) - 
2f C(5) - SC(3) 


1.022090029 


m = 8 


4£o+5C(9)-^C(3)-£C(5)- 
C(5)C(3) - C(7) (l + £ ) 


1.010557246 
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m = 9 


34020 9450 C(7)C(3) 

5C(9) - £ C(7) - ^C(3) 


K(5) 2 + 

- §k(5) 


1.005133570 


m = 10 


3iSo+6C(H)- 

K(5) 2 - 

c(9) (i + f ) 


9450 ^( 3 ) 945 

(7) - C(7)C(3) - 


1.002522063 


5(3, m) 


Power of n 


Closed form 


Approx. value 


m = 4 


S - ? " ?C(3) - K(3) 


+ 3C(5) 


1.628620203 


m = 5 


iio^fi4 _ TC(3)-4C(3)- 
1C(5) 


K(3) 2 + 


1.180103635 


m = 6 


3S-T^+ 4 C(7)-^C(5)- 
IC(3) 2 + lC(5)-fJC(3) 


£c(3)- 


1.072362484 


m = 7 


7T 8 7T b 7T 4 A/Q\ 

4200 2520 60 ^\ ' 

C(5)C(3) - C(5) + f ) 


K(3) 2 - 

+ 6C(7) 


1.032351029 


m = 8 


2800 1260 

C(5) (S + £ 

c(7) (! - * ) 


C(3) (fg) + A) 
) - |C(5)C(3) 
+ 5C(9) 


+ 


1.015179175 


5(4, ra) 


Power of n 


Closed form 


Approx. value 


m = 5 


7T b 7T 4 ll7T^ >-/Q\ 117T 2 Q"\ 

540 810 216 36 

iC(3) 2 + ±R(5) 


1.310990854 


m = 6 


llTT 6 7T 4 A/o\ / 7T 4 , l7T 2 

3240 80 >V J ^90 "r 6 

+i|C(3) 2 + C(5) 


T 36 J T 

)+4C(7) 


1.103348021 


m = 7 


7T S llTT 4 7T b Vf^ 2 

4200 3240 ' 2430 2^°/ 

C(3) (^-i) - C(3)C(5) 

c(5) (§-¥) + f C(7) 




1.043816710 



10 
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" C(3) (I + g) - 



ll7T M _ 577° 

25200 4536 



^C(3) 2 - §C(3)C(5) 
C(5) + T + §) 
C(7) (j§-^)+5C(9) 



+ 



1.020093103 
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